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ABSTRACT 


In the present work optimization of reinforced 
concrete elastic grid roofs is carried out. The objective 
function chosen is the total cost of materials with const- 
raints on maxiraum shear stress and maximum compression 
reinforcement in the sections. The problem is formulated 
as a mathematical programming problem and is solved using 
sequential unconstrained minimization technique with interior 
penalty function approach. Tariable metric algorithm is 
used for unconstrained minimization. 

The grid structure is analysed by stiffness method 
of structural analysis and the resulting equations are solved 
by Choleski's decomposition method. 

Parametric study is made by varying the nimiber of 
beams in the grid. An empirical formula for obtaining the 
optimum number of beams is suggested. 



CHAPTBE 1 


lETROIiUCTIOl 


1 . 1 G-EEERAL 

Concrete and steel are widely used materials in 
civil engineering construction works . With the advent of 
digital computers and development of optimization techniques 
substantial saving of these construction materials can be 
achieved. Present work attempts to obtain the minimum cost 
design of large elastic reinforced concrete grid system via 
mathematical programming. 

designers are sometimes called upon to provide 
designs to cover up large floor areas with the following 
conditions . 

1 ) The entire floor area to be clear of any 
structural obstruction 

2) The upper floor to be flat and 

5) The head-room to be maximum. 

The first condition can be achieved by eliminating 
intermediate columns and providing deep main beams- This 
solution does not fall in line with the third condition but 
can be improved upon by the use of inverted beams. This 
again does not satisfy the second condition. 
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The solution to this problem is the use of grid 
system, constructed of reinforced cement concrete beams with 
thin slab on top of it. 

In a grid system there are usually two sets of beams 
running in two different directions and are rigidly joined at 
the points of intersection. In general there can be more 
than two sets, their respective elements need not be parallel 
and the angles of intersection are arbitrary. However, this 
general form of grid is rarely encountered in practice. The 
grid systems normally used are orthogrid and diagrid which 
are shown in Figures 1.1 and 1.2 respectively. In orthogrid 
there are two sets of beams running parallel to the supports 
and which intersect each other at right angles. In diagrid 
system the parallel beams running in two directions are 
inclined to the spans. The angles of inclination of the two 
sets with the spans need not be same . 

The grid system is a statically indeterminate 
structure and the degree of redundancy depends upon the 
niunber of beams used in each direction. The load on the grid 
is essentially the dead load and the live load. The load 
shared by a beam will depend upon the spacing of the beams . 
Thus the spacing of the beams becomes an important parameter 
from the point of view of the designer who is always 
conscious to provide a minimum cost design. Keeping in mind 
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that the grid system is highly indeterminate and the objectiv 
of the designer to have the minimum cost design, this problem 
is most suitable to be formulated as optimum design problem, 
the solution of which can be tracked down with the help of 
digital computer. 

1 .2 STATBMEITT OF THE PEOBIEII 

Present investigation is undertaken to find optimum 
elastic design of reinforced concrete grid roofs. The para- 
metric study for the mmnber of beams leading to the minimum 
cost design has been made. A grid square in plan with two 
sets of equally spaced beams , having same spacing in the two i 

directions and running parallel to the supports is considered . : 
The grid is considered to be fixed on supports . The width 
and depth of the beam sections in both directions is kept same J 
The slab is considered to be simply supported over the grid. 

The grid is analysed with stiffness method as explained in 
Chapter 2. The cross sections are designed according to the j 
elastic design philosophy which is explained in Chapter 3 . The 
formulation of the mathematical programming problem is explain^ 
in Chapter 4. Chapter 5 presents results, discussion and I 

conclusions. i 
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1 .3 EEVIEV/ OP RELATED TOPICS 
1,3.1 Analysis of Grids 

lazarides (21) used deflection compatibility equations 
at the intersection of beams to analyse a grid system. This 
turned out to be a tedious and impracticable method. Moment 
distribution method was applied for grid analysis by Halasz 
(13) and also by Ewell, Okubo and Abrams (4). This method also 
inTol'T'ed lengthy computations. Guyon (12) obtained a solution 
by representing the grid system with an anisotropic plate 
where, the members were assumed to have negligible torsional 
resistance. This solution was generalized by Massonnet (1,23) 
to take the torsional rigidity into consideration. This 
equivalent plate method was applied by Taraporewalla (30) for 
designing grid systems. Hendry and Jaeger (14) proposed to 
replace the cross-beams of the grid system by a uniformly 
distributed medium. Here torsional rigidity of beams is taken 
into account but that of the cross-beams is neglected. The 
analysis is done by solving the differential equation of 
loading by harmonic analysis, lightfoot and Sawko (22) applied 
generalized slope-deflection method for the analysis of grid 
system with the help of computer. Reilly' (27) incorporated 
the effect of warping in the analysis of grids by stiffness 
method. Work in' plastic analysis of grids has also been done 
by Hodge (17), Heyman (15,16) and Grigorian (9,10). 
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In the present work stiffness method of structural 
analysis is used to analyse the elastic grid. This method is 
amenable to handle large redundancy and is simple to programme 
on computer. 

1.3.2 Optimization of Grids 

Optimization of simply supported plane steel grill- 
ages carrying uniformly distributed loads was studied by 
Grundy (11). Here weight function was obtained from mid-span 
bending moment by limiting extreme fibre stresses or iiltimate 
moment or maximum deflection to a permissible value. Assuming 
the weight function to be differentiable, optimum values of 
spanswere found out by partial differentiation. 

Optimization of a simply supported reinforced 
concrete square grid roof was studied by Harayanan (24). The 
analysis was carried out by representing the grid with an 
equivalent orthotropic plate . The section was designed by 
plastic method for mid-span bending moment, for the optimal 
design, section of beam, depth of slab, percentage of steel in 
beam and slab and spacing of beams were found out. 

Ramesh and Karve (26), carried out parametric study 
for orthogonal grids with different boundary conditions and 
with varying number of beams . The system was analysed as a 
stiffened plate adopting a conforming rectangular 
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plate element and using finite element technique. The para- 
meters studied are, the ratio of span to plate thickness and 
the ratio of span to depth of beams . Ratio of span to width 
of beams was taken equal to 24 in all cases. Charts are 
presented for different grids of given aspect ratio and given 
number of beams running in each direction, for a given grid 
the values of the above mentioned parameters, which achieve 
given effective thickness with minimum volume of material can 
be found out from the charts. 

Keshav Rao (19) used fully stressed design method 
for finding an optimal layout for grid. In this method 
starting with a unifoim depth of slab the depths at various 
nodes of the slab are changed continually to obtain fully 
stressed design in each cycle. After a certain number of 
iterations the change in nodal depths from one cycle to next 
cycle becomes negligible. At this stage the ridges obtained 
by joining the nodes of large depth are interpreted as beams. 
From this a practical layout of beam system is obtained 
showing substantial saving of material. 

Kwlecinski and Kleiber (20) studied optimization of 
simply supported rectangular grids composed of two recti- 
linear families of perfectly plastic bars subjected to 
transverse load. Grid bars were assumed to be capable of 
resisting torque. Optimum configuration was found based on 
minimizing weight of the structure . 
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The present work studies the optimum spacing of 
beams leading to the minimum cost design of the reinforced 
cement concrete grid system. 

1 . 3.3 Optimization Techniques 

Optimization techniques for constrained problems 
can be broadly classified into two groups. One, direct methods 
and other, indirect methods. Of the direct methods, 
Zoutendijk's method (31) of feasible directions is used for 
problems with non-linear inequality constraints. The gradient 
projection method by Rosen (28,29) is well used for problems 
with linear constraints but is ineffective for non-linear 
constraints. The indirect method of penalty function (8) is 
simple to apply for the constrained problems. There are two 
approaches for the penalty function method. One is exterior 
penalty function method and other is interior penalty function 
method. The latter has certain advantages over the exterior 
penalty function method. In the approach by interior penalty 
function, 1) a feasible design point is always at hand, 2) 
improved designs are obtained in a sequential process, and 
3) the constraints become critical only at the optimum. Hence 
in the present work interior , penalty function approach has 
been used. This method requires the use of an unconstrained 
minimization technique. 
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Steepest descent method was used in early days for 
unconstrained minimization ^ But for eccentric functions this 
technique faces difficulties in convergence. Bletcher and 
Reeves (7) developed the conjugate gradient technique popularly 
known as Rletcher-Reeves method. Powell (25) introduced the 
conjugate direction method in 1964 for unconstrained minimi- 
zation. Devidon (3) invented the variable metric method to 
carry out unconstrained optimization. This was further 
sharpened by Pletcher and Powell (6). The method is also known 
as DPP method. This technique is powerful and has been 
successfully used for optimizing highly distorted functions. 

For present work Sequential Unconstrained Minimization 
Technique (SUMT) of Piacco and McCormick (5) is used which is 
an interior penalty function method. The PPP algorithm is 
used as the technique for unconstrained minimization and 
linear minimization is carried out by polynomial fitting. 
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CHAPTER 2 
AMIYSIS OP GRED 


2.1 IITRODUCTIOF 

Stiffness method of structural analysis is used for 
analysing the grid in the present work. Each beam intersection 
is considered as a node with three degrees of freedom viz., 
the transverse deformation and the rotations about two 
directions. Advantage of the symmetry of structure is taken 
to reduce the size of the problem. To seek the solution of the 
stiffness formulation Choleski's decomposition technique is 
used. 

The materials considered to be used for the const- 
motion of the grid are; 1) Concrete of Ml 50 grade and 2) The 
reinforcement to be of Grade I mild steel. The properties of 
concrete and steel are in accordance with IS-456-1964. 

2.2 ASSUMPTIONS 

Hollowing assumptions are made in the analysis of 
grid for the present work. 

l) The intersection of the beams i.e. the nodes are 


rigid joints. 
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2) Eeinforced cement concrete is homogeneous, 
isotropic and linearly elastic material. 

3) The stiffness of iDeam is obtained from the 
overall section of the beam neglecting the steel. 

4) The value of the coefficient P in the expression 
for the torsional stiffness (pG-b^T) is kept equal to 0.3 for 
all values of b/D. 

5) Poisson’s ratio for concrete is assumed equal to 


zero. 


6) Axial deformations in the beams are neglected. 


2.3 C0ESIDEEA.TI0N OP SYMMETRY 

Only a quadrant of the structure need be analysed by 
taking the advantage of symmetry about two directions. While 
solving the simultaneous equations, the condition is imposed 
that the slopes of beams crossing the line of symmetry are 
zero at the line of symmetry. 

figures 2.1 and 2.3 show a quarter of the symmetric 
grid structure when the number of beams running in each 
directions (Nb) is odd and even respectively. In the case 
when ITb is odd the lines of symmetry pass along the central 
beams in two directions. In such case symmetry is taken into 
account by considering half the stiffnesses of the beams 
running along the lines of ssnmnetry. When Fb is even 




"•H '■ “‘v ' - 


















Half span 


FlG-2'3 GRID WITH EVEN NUMBER OF BEAMS IN EACH 
DIRECTION ■ 

^ . ' ■■ ; ./ ' ■; O, ,/, ■ TO"" ir. i‘„ , 

^.Load in this area to be 

shared by this beam' 


Load in this area to be 
shared by this beam 


Beams supporting 
the two way slab 


distribution of load on beams 
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half the length of elements crossing the line of ssnmnetry 
comes into picture . Here the points where the elements are 
cut are taken as nodes and the stiffnesses of these elements 
are calculated according to their length effective in the 
quarter . 

2.4 LOiD ON BEAMS 

As the spacing of the beams is same in both the 
directions the amount of load transferred by the slab, which 
acts as a two way continuous slab , to each element is same 
(see Eigure 2.4). This load consists of live load on the 
slab and the dead weight of the slab and is taken to be 
uniformly distributed over the elements. This load plus the 
self weight of the beams constitutes the total load on beams 


2.5 EORCE DEEOEMATIOH EEIATIOKTS 

The grid is subjected to transverse load only and as 
previously stated axial deformations of the elements are 
neglected. Hence the deformations which need be considered at 
each node of the element are , 

the vertical deflection = u 

the flexural rotation = 0 

and the torsional rotation = 0 
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Corresponding forces at each node of the element are, 
the shear force = V 

the "bending moment = W 

and the torsional moment = T 

Element notations for the forces and deformations 
of near and far end of the elements are shown in Figure 2.5* 

The superscripts If and E denote the near and far end quantities 
respectively. The near and far ends of the elements running 
in X and Y directions are shown in Figure 2.6. The glo"bal 
notations for forces and deformations are same as the element 
notations for elements running in X-direction. 

let us denote the nodes and elements as follows: 

"tin 

(i,j) is the node on i row and j column of the 

nodes . 

x(i,j) is the element running in X-direction and 
having it's near end at node (i,j). 

y(i,j) is the element running in T-direction and 
having it's near end at node (i,j). 

where i = 1 , 2 , 3 » • • • 

Fow the final forces at each node of elements x(i, j ) 
and y(i,j) can "be written as (refer Figure 2.6). 

„fF ^ , F „ ^ F „ 

^x(i,j) ^x(i,3) ^(i,j) ^ ^^x(i,o) ^(i,j+1) 


f 

x(i, j) 


( 2 . 1 ) 




FlG-2-5 ELEMENT AND GLOBAL NOTATIONS FOR FORCES AND 
DEFORMATiONS 




Element ! Element 


AND FAR ENDS OF THE ELEMEN 
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In the above equations the suffix x(i,j)5 y(ijj) 
denote the quantities corresponding to these elements. The 
suffix (ijj) denotes quantities corresponding to node (i,j). 
The superscripts If and P denote quantities corresponding to 
near and far ends of the elements respectively. Other 
notations for an element are as follows. 

f is the final force vector in global notations 

given by 
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f is the fixed end force vector in element 
notations . 

k is the direct stiffness matrix of the element. 

For the near and far nodes the direct stiffness matrices are 
given hy . 



1 

H • 

6EI 

l2 

0 


4EI 

1 

6EI 

■ 

0 

III 

6EI 

1^ 

12BI 

y 

0 

■p 

and k s 

6EI 

1^ 

12EI 

G 

0 


0 

0 

1 _ 


0 

0 

GJ 

1 


Here 

hl^ 

I = is the moment of inertia of the element section 

b is the width of the section 
D is the depth of the section 

1 is the length of the element i.e. the spacing of beams 
G-J = GPb D is the torsional, rigidity of the section 
E is the Young's modulus of elasticity for concrete 
G is the modulus of rigidity for concrete 

ck is the carryover stiffness matrix of the element. 
The carryover stiffness matrices for near and far node defor- 
mations are given by 
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ck- 




and 


1 E" 
ck 


2EI 

1 


6EI 


2EI 

1 


6EI 
2 




0 


6EI 

l2 

12EI 

0 


6EI 

l2 

12EI 

0 


0 

0 

G£ 

1 

0 

0 

1 


U is the deformation Tec tor in global notations 


given by 


9 

u 

0 


t is the transformation matrix, which by premulti- 
plying to a deformation vector U/ . converts it into element 
notations for elements y(i,j). 

t = (t) = transformation matrix which by 

p remalt iplying to a force vector of element y(i,j) in element 
notations, converts it into global notations. 
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t and t are given by 



0 

0 

1 

- 

0 

0 

-1 

t = 

0 

1 

0 

III 

0 

1 

0 


-1 

0 

0 j 


1 

0 

0 


2.6 EQUATIONS OF ' EQUILIBRIIM 


The equations of equilibrium of the forces at the 
node (i,j) can be written as 


^(i,j) ^(i-1,o) ^y(i,3) 


W/ - -X 


(2.5) 


where W/. is the external load vector at node (i,o) in 
j 3 >' 

global notations. 

Using equations (2.1) through (2.4) » equation (2.5) 
can be written as 


N * N 

°^x(i,j-1) ^(i,3-l) ^ ^ ^ ^(i-ljj) 


+ 


k: 


F 


* F * N ,N 

/- .^\4"tk/'.^ -xt + tk/- .\t + k/ . 

x(i,3-l) y(i-1,3) t(i» 3) x(i,j) 


U/. .X 

(i»j) 


*_,F 


^ ^ “^(1,3) ^ hl+1,3) * “(1,3+1) 

* 

(i,j) “ ^x(i,j-l) “ ^ ^y(i-1,j) 


P 




=,fN 


w 4,fF +*.fP ■ +*^f]T _ 

- W/_. .X - f./. . .x t f„^,. ^ - t ^x(l,2) 


( 2 . 6 ) 
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again this is written as 


=*"x(i,3-1) ’^(1,3-1) °‘^y(i-1,3) '^(1-1,3) + ^(1,3) '^(1,3) 


+ CkJ,. U,_. . . + Ckl/. U. 


_F 


y(i,3) (i+ljj) x(i,3) (i,j+1) 


= ¥ 


(ij 3) 


(2.7) 


where , 

modified carryover stiffness matrix of element y(i-1,3) for 
it's near end deformations is given by 


CK- 


-U 


* N 

^ "Vi- 1 , 3 ) 


y(i-1,D) 

modified carryover stiffness matrix of element y(i^j) for 
it's far end deformations is given by 


GK^ 


* m 


y(i,3) 


= t ck; 


7(i,3) 


combined direct stiffness matrix at node (i,j) is given by 

1* * I* . * F , F 

Ki. .X = k /. . + t k / . , .X t + t k /. .\ t + k^/. .n 

(i,y) x(i,j-l) y(i-1,j) y(i,j) 2 :(i,o) 

modified external load vector at node (1,3) is given by 
^(i, 0 ) ^(i,3) ^x(i,5-l) ^y(i-1,3) ^y(i,3) 


„fF 

x(i,5) 
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Equation (2. 7) is written for each node and all the 
equations are combined and written in the form 

AX = B (2.8) 

where A = the stiffness matrix or coefficient matrix of 
the system 

X = the deformation rector of the system 
B = the modified external load vector of the system. 
Equation (2.8) is the set of simultaneous equations 
with the deformation vector X as unknown. 

2.7 SOLTJTIOlSr OE THE SimiTAHEOHS EQUATIONS 

The actual numbering of nodes is done as shown in 
Figures 2.1 and 2.3. This pattern gives a banded stiffness 
matrix of the system. The structixre of the coefficient 
matrix A as given by equation (2.8) is shown in Figure 2.7. 
The solution of the equation (2.8) is obtained by Choleski's 
decomposition which is explained below. 

The matrix A is sjnmmetric and can be written as 

A = S^S (2.9) 

where S is an upper triangular matrix. 

Using equation (2.9) in equation (2.8) we get 

T — 

S SX = 


B 


( 2 . 10 ) 




Half band -width 


Full band width 
P-2p-1 


Diagonal 

elements 


FiG-2-7 STRUCTURE OF THE COEFFICIENT MATRIX A 


Band width 






Dt agonal" 
elements 




OF THE DECOMPOSED MATRIX S 
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putting SX = Y in atove equation we get, 


( 2 . 11 ) 


rp— _ 

S Y = B 


( 2 . 12 ) 


Once the matrix: S is obtained, equation (2.12) can 
be solved for Y by forward substitution and using Y in 
equation (2.11), X can be found out by backward substitution. 

Mow to find out the matrix S let us consider the 
elements in the upper triangle of A, We denote, 

a. . = element of matrix A lying in i"*^^ row and 
column 

s . . = element of matrix S lying in i"^^ row and 

column 

!3? "fell. "t/ii 

s! . = element of matrix S lying in i row and j 
1 s 0 

column. 

Brom equation (2.9) the elements of matrix A can be 
written in terms of the elements of S as, 


®i,1®1,5 ®i,2®2,j ••• ®l,l®i,3 


i j 


(2.13) 


= S-, -s. .+ s^ .s^ ,+ ... + s. .s. . 


1 IS 3 


(2.14) 


4,3 = 
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For diagonal elements putting j = i in equation 
( 2 . 14 ) we get, 


^i,i 


2 , 2 
,i ®2,i * 


. + s. . 


( 2 . 15 ) 




From equations (2.14) and (2.15) 


s. . 


®i-1 ,i®i-1 , 1 


ana = (=-1,1 - =1,i - 4,i - ••• - 


(2.16) 


(2.17) 


Using equations (2.16) and (2.17) elements of S 
can be found out . 

The structure of the matrix S is also banded with 

band v/idth equal to half the band width of matrix A. This 

can be easily explained as follows . Consider an off -band 

element a. . of matrix A in the first row 
' > J 


••• ° 


J > P 


!2.18) 


from equation (2.16) we have 


j > P 


from equation (2.18). 
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Extending the argument for elements in other rows 
it can he proved that 


s-i i = 0 D - i 2: P 

\ 

The structure of matrix S is shown in Figure 2.8. 


2.8 FORCES IF THE ELEMENTS 

After the equation (2.8) are solved for deformations 
the final forces at each node of the elements can he found 
out with the help of equations (2.1) to (2.4). 
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CHAPTER 3 

PESIGH OP SECTIOHS 

3.1 IFTRODUGTIOH 

Working stress design method is used for designing 
the sections of each element of the grid in the present work. 
Forces at both nodes and at the centre of the element being 
known from analysis, the maximum force is taken for designing 
the element section. The section is designed for three forces 
viz . 5 the bending moment , the torsional moment and the shear 
force. Specifications from IS code of practice (18) are 
followed for the design. 

3.2 ISSHlVlPTIOSrS 

Following assumptions are made in designing the 

section. 

1 ) Plane cross-sections before bending remain plane 
after bending. 

2) For concrete stress is directly proportional to 
strain i.e. Young’s modulus of elasticity is constant at all 
stresses. 

3) Concrete area on tension side is ineffective. 

4) Reinforcing steel is concentrated at the 


centroid of the steel area. 
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3.3 DESIGN QP silB 


The slah mins over a mimber of beams running at 
right angles to each other. Hence it acts as a two way 
continuous slab and is designed accordingly. The design is 
done according to "The method 5" as described in IS-.456-.i964. 
In this method tables are presented which give coefficient* 
for maxiiaum bending moment in a two -way continuous slah with 
given aspect ratio. 

Minimum thickness of elab Is taken as larger of the 
following: 


1 ) 


scan 


10 cms 


2 ) = to 

where t is miniTnum thickness of the slab, 
smin 

It is' seen that for a given spacing, specified 
minimum thickness of slab is sufficient to carry the load 
safe>.y. This is true at all the spacings considered in the 
present work. 

Mi ni muTtv percentage of reinforcement in the slab 
is taken to be equal to 0.155^ of the cross-section. 
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3,4 DESIGN OE BEAM EOH PIEXIIRE 

The heam is designed for flexure such that the 
section is either singly or doubly reinforced. When the 
lending moment at the section is less than or equal to the 
moment of resistance of the balanced section (MRb) the 
section is designed as a singly reinforced beam and in other 
case i.e. when the bending moment is greater than MRb the 
section is designed as a doubly reinforced beam, following 
sections explain briefly the design of singly and doubly 
reinforced beams. 

3.4.1 Design of a Singly Reinforced Section 

Referring to figure 3.1 the position of neutral 
axis and moment of resistance of a singly reinforced section 
can be obtained from equations 


bn 

2 

= mA (d 
s 

- n) 

(3.1) 

MR 

s s 

- -) 

y 

(3.2) 


where 

n is the distance of neutral axis from the extreme 
fibre in compression 

d is the effective depth of beam 
m is the modular ratio for steel and concrete 








Beam section 


Stress diagram 


FIG' 3-1 SINGLY REINFORCED SECTION 


Neutral 

axis 


Stress diagram 


Beam section 




DOUBLY REINFORCED SECTION 


Neut ral 
axis 




MR is the moment of resistance of the section 
A is the area of steel in tension 
and al is the maximum permissible stress in steel. 

O 

To find the steel area A such that MR is equal to 

s 

the external bending moment M, is a trial and error process. 
These trials are eliminated in the present work with the 
following approximation. 

A , X M 

A _ 

s ~ MRb 

where 

A V is the area of tension steel for the balanced 
sb 

section. 

The steel area obtained from the above approximation 
is slightly more than the actual required steel. Bht the 
difference is very small and can be neglected for practical 
purposes. 

The minimum percentage of steel in the beam is kept 
equal to 0*3% of the overall section. 

3 . 4.2 Design of a Doubly Reinforced Section 

Referring to Figure 3.2 and taking the moment of 
compression force about the neutral axis we get, 
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i(d 


(m - 1 )A cr ' (- 
sc c ^ 


-)(d - c) 


= MRb + M. 


where 


a' n 

-I (d - -I 


(3.2) 


is the moment of resistance due to concrete area, 


c 


(3.3) 


is the moment of resistance due to compression steel, 

n is the distance of the neutral axis from the 
c 

extreme compression fibre in a balanced section, 

cr' is the maxinrum permissible bending compressive 
stress in concrete, 

is the area of compression steel, 
and c is the cover for compression steel. 

Equating Iffi. and M we have 


MR = K = MEb + 1, 


(3.4) 


BM - MRb 


(3.5: 


MRb can be computed for a given section by using 
equation (3.2). Furthermore, using equations (3.3) and (3*5 


can be found out. 



53 


How equating moment of areas about the neutral 
axis we get, 

~ n^) = (m - - c) + bn^ (5.6) 

where from the only unknown A i.e. the tension steel area 

s 

can be computed . 

The value of A is subject to a maximum of 4% of 

o C 

the gross section as per IS code of practice. 


3.5 IBSIGU BOR SHEAR AND TORSION 

The shear stresses in a section due to the shear 
force and torsional moment are given by 




JL 

h a 


and 

where 


12 = 


T(3 + 
b^D 


is the shear stress due to shear force Y 
is shear stress due to torsional moment T 

a is the lever arm of the section. 

Now, if q-| + 12 <I l’ ■fcb.e maximum permissible shear 
stress in concrete, nominal shear steel is provided. 
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if + q 2 4(q') the section is discarded and 
if q' -c, q^ + q^ -<f 4(q') shear steel is designed as follows 

For shear force 


q^b 


ssl 


cr 


ss 


and for torsional shear 


T 


^ss2 


0.8 o'X.Y. 
s 1 1 


where -^gg-j 4s the steel area per unit length in all legs 
of stirrups for shear force 
Ags 2 steel area per unit length in all legs 

of stirrups for torsional shear 


and 


o' is the permissible tension in shear steel 
ss ^ 

and are the sides of the rectangular stirrups. 


In addition to stirrups longitudinal steel is 
provided such that it has the same volume per unit length 
as the volume contained in all legs of stirrups resisting 
tors ional shear. 
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CHAPTER 4 

EOEMUIiATIOH OP THE MATHEMATICAL PROGRAMMING PROBLEM 


4.1 INTRODUCTION 

Once the method of analysis and the design 
philosophy are decided upon, the optimum design pro-hlem can 
he cast as a mathematical programming problem which has the 
following form. 

Minimize P(D) , D = 

% 

Subject to conditions 




i = 1, 2, ...» NC 


where 


and 


PCD) 

N 

gjD) 

NC 


is the objective function 

is the number of design variables 
ijli 

is the i constraint 

is the total number of constraints. 
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The details regarding the design variahles, 
ohjectiTe function, constraints and the technique for opti- 
mization are exjplained in the folio-wing sections. 

4.2 THE DESIGN VAEIABDES 

Design variables in an optimization problem are the 
independent quantities for which values are to be assigned 
in producing a design. For designing a grid roof of given 
size and with a given niunber of beams the quantities which 
are required to be considered independently are, the slab 
thickness, width of the beams and depth of the beams. In 
the present work the slab thickness is always kept to 
minimum as required by code in order to keep the dead weight 
on the grid at it's minimum. Hence only two quantities viz. 
width of the beams and depth of the beams form the design 
variable vector. 

4.3 THE OBJECTIVE EIJNCTION 

For a grid with given number of beams i.e. for a 
given spacing of beams the cost of slab is constant as it's 
thickness is kept constant. Hence only cost of beams is 
going to change as the design variables change. So the 
objective function considered for optimization is the cost 
of beams, the expression for which is as follows; 



37 


P(b,D) = m>(Kb' + 1) -b D(21 - b)C 

c 

""si ^soi + (^ssli * 2^ss2i’^r^l) 

where F is the cost of beams in the grid 

cost of concrete per unit volume 
Cg is the cost of steel per unit volume 
Ml is the total number of elements in the grid 
= Fb(m) + 1)2. 

The cost of materials are taken from Delhi Schedule 
of Rates (2). Cost of concrete is taken as Rs. 166.6 per m^. 
It does not include the cost of formwork. The cost of 
reinforcement is taken as Rs. 2.2 per Kg which includes the 
cost of bending, binding and placing in position. 

4.4 THE COKSTEAIKTS 

Constraints in working out a design are the 
conditions which when violated give an infeasible design. 

In a design problem certain limitations are imposed on the 
behaviour of the structure to get a safe design. These 
limitations are called as the behaviour constraints. Usually 
there are still other limitations required to be considered 


+ C 1 
s 


E 

i=1 
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so as "to olDtain a practicable design. These are termed as 
side constraints. For the present work following constraints 
are considered to get a safe and practicable optimal design. 
The behaviour constraints are % 

( 1 ) limitation on maximum shear stress 

= ^l(j) ^ ^2(3) - 20 < 0 (4.1) 

3=1,2,..., NED 

I 

(2) limitation on the maximum compression steel area 


^ 2 ( 3 ) 


■^ 60 ( 3 ) 



(4.2) 


3 = r, 2, ...» NED 


where NED is the number of elements designed. 

The side constraints are as follows; 

(3) limitation on the minimum width of beam. This is taken 
as the critical of the following; 

a) limitation on the absolute minimum width of beam 


g^ = 15 - b < 0 (4.3) 

b) limitation on the width of beam so as to avoid 
a slender beam 

= 1 - 30b 0 (4.4) 

(4) limitation on the maximum width of beam 


§4 


= b - 


0 


(4.5) 
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(5) limx'ta'tion on "ttie non— nega'tivi'fcy of th© depfli 

= - D < 0 (4.6) 

4.5 TECHNIQUE FOR OPTIMIZATION 

The constrained minimize.tion problem is transformed 
into a sequential unconstrained minimization problem with 
the use of interior penalty function. The penalty function 
used for present work is 

r) = P(E) - X t — (4.7) 

i=1 g^CD) 

where jZ!(D, r) is the penalty function 
and r is the penalty parameter. 

Sequential unconstrained minimization comprises 
of the following steps t 

1 ) With a feasible design point and an initial 
value of r minimize 0iDy r) to obtain 

2) Check the convergence for to.be optimum 
design. If satisfied terminate otherwise go to next step. 

3) Reduce the value of r by r hr where h is 
the reduction factor. 

4) Compute a new starting point and repeat the 
minimization process from step 1 , 
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To obtain a feasible starting point in order to 
initiate the optimization algorithm, only a few trials were 
required in the present problem. Initial value of r was 
chosen such that at the starting point 1^ the penalty 
imposed is equal to the objective function value. In the 
sequential process of minimization r was reduced by the factor 
h = 0.1 till it reached a value equal to 0.01 and then the 
algorithm was terminated . The starting point chosen for 
the reduced r was the minimum point obtained for the previous 
value of r , 

The unconstrained minimization of 0(I>, r) was 
carried out by the variable metric (IFP) method. Brief 
outline of this method is presented below. 

1 ) Set q = 0. 

2) Start with an initial design point and an 
initial positive definite matrix which is taken as 
identity matrix in the present work and set the initial 
direction vector 

So = - HoV^(\. 

3) Compute \ “A 

where is the minimizing step length in the direction 
giving 


as the point at minimum 
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4) Check; the coayergence for 4o he the 


optimiim design point < 
to next step* 

5) ^.Compute 


if Satisfied terminate otherwise 


* s ■ 

H = H + a - 

q+1 q q 


■ (H Z )(H-Z ) 


T 


Z^H z 
q q q 


where Z^ = y , r) -.\70(D^, r) 

6) Compute = "Vi Vl » 

7) Set q ^ q+1 and repeat from step 

ihe convergence criterion used in step 4 is that 
the percentage reduction in the penalty function value 
during two successive linear minimizations is less than 1?S. 
The gradient evaluation is done by using the forward 
difference technique and for linear minimization quadratic 
interpolation is used. The quadratic was fitted to an 
accuracy of 15?. 
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CHAPTER 5 

RESULTS, LISGUSSIOH AED CONCLUSIONS 

5 . 1 GENERAL 

The application of the formulation described in 
the preTious chapter is explained with three illustrative 
examples in the following sections. The results presented 
in this chapter are obtained on IBM 7044/1401 system at 
Indian Institute of Technology, Kanpur. 

The numerical values pertaining to material 
properties used for solving illustrative problems are as 
follows ; 

1 ) Grade of concrete Ml 50 

2) Reinforcing steel M S Grade I 

3) Density of reinforced cement concrete 2400 Kg/m^ 

rz 

4) Density of reinforcing steel 7880 Kg/w 

5) Maximum permissible bending compressive 

2 

stress in concrete 50 Kg/cm 

6) Maximum permissible tensile stress 

2 

in reinforcement 1400 Kg/cm 

7) Modular ratio for steel and concrete 18 
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5.2 ILLUSTEAIIYE EXAMPLE 1 

In this example a grid with overall size 
18.82 M X 18.82 M (Grid 1 ) is considered to find the optimum 
spacing of beam leading to the minimum cost of the roof. 

For illustration let us consider a case when the number of 
beams running in each direction ai'e even, say equal to 8. 

The quarter of this structure (Figure 2.3) has 24 nodes. 

Thus the total number of degrees of freedom turn out to be 
72. The spacing of beams is 209.11 cms centre to centre. 
Therefore, the slab thickness comes out to be 10 cms i.e. 
the Biinimum to be provided. This gives a dead weight of 
240 Kg/m on the grid. In addition to this dead weight, a 
live load of 75 Kg/m on the slab and the dead weight of 
the grid (which varies as the design changes) are also 
considered. 

As indicated in Chapter 4, the optimum design 
problem is a two variable problem with the width and depth 
of beams as design variables. Therefore L = (b, D) . 

The width of the beam is allowed to vary in the 
range of values such that the beam does not become a slender 
beam on one hand and the beams still maintain a spacing: 
between them. In other words, for the example considered 

209.11 _ _ 209.11 

50 ^ ^ 2 


6.97 


104.5 
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In order to make the lower and the upper hounds 
worked out above as practicable the constraint actually 
imposed on the width is 

15 b :< 60 

The depth of the beam is only constrained to be 
non-negative . 

The difference in the levels of tension reinforce- 
ment in X and Y direction beanB is taken into account by 
providing a cover of 4 cms and 6,5 cms in the two directions 

respectively. The cover for compression reinforcement is 

*t/ln 

assumed to be 1/10 of the effective depth of the section. 

The cross-section of the grid between two nodes is 
designed for maximum bending moment, shear force and 
torsion and reinforcing steel is accordingly provided. 

Hence for each design twenty cross-sections are to be 
detailed in each direction. For each cross-section the 
maximum compression steel and maximum shear stress are 
constrained. Thus the optimum design problem turns out to 
be having a total of 83 constraints. 

The problem has been solved with a starting point 
(45.0, 130.0) having a total cost of Rs. 62,007. The 
proposed optimum design comes out as (15.02, 129.39) having 
a total cost of Rs. 26,457. The computer time taken for 
this problem is 349 secs. 
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A parametric study with different number of beams 
ranging from 2 in increments of one upto 10 beams in each 
direction has been made. The results of the parametric 
study are shown in Table 1 . A glance on the table shows 
that optimum number of beams turns out to be 4 with a spacing 
of 376.4 cms. The total cost of bhe optimum design is 
Rs. 20 , 763 . A plot of number of beams in each direction 
■versus total cost of roof is shown in figure 5.1. 

As a check for the optimuni design to be local or 
global minimum, minimization is carried out with one more 
starting point (35.0, 100.0) which gi'vesthe optimum cost of 
roof as Rs. 20,766 at a section (15.05j 124.37). Since in 
the two cases the optimum design reached is similar and so 
also the optimum cost , it can be said with some degree of 
confidence that the proposed optimum design is the global 
minimum of the problem. 

5.3 ILIUSTRATIVE EXAMPIE 2 

A second grid with o'verall size 30 II x 30 M (Grid 2) 
is optimized for the purpose of studying the optimum spacing 
of the beams . for illustration consider the grid when the 
number of beams are odd, say equal to '5. The quarter of 
this structxxre has 9 nodes resulting in 27 degrees of 
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Table 1 



Grid 1 

; Size 18.82 

1 X 18.82 

M 


o 

« 

, o 
i ^ 

1 

1 Starting design point 

; Proposed 

optimum design point 

beams in each 
direction Ub 

I 

5 b cms 

I 

j 

i 



D cms 

Cost of 
roof 
(Rs.) 

! ! 

1 b cms J 

1 f 

i 1 

1 ! 

1 

f 

! 

I) cms J 

f 

t 

t 

Cost of 
roof 
(Rs. ) 

10 

40.00 

150.00 

62,476 

15.08 

145.43 

31 ,021 

9 

40.00 

120.00 

52,789 

15.00 

122.13 

27,136 

8 

45.00 

130.00 

62,007 

15.02 

129.39 

26,457 

7 

45.00 

135.00 

57,187 

15.04 

133.00 

25,034 

6 

50.00 

160.00 

55,874 

15.00 

113.69 

22,825 

5 

50.00 

150.00 

54,311 

15.00 

140.90 

22,105 


55.00 

150.00 

52,178 

15.08 

130.85 

20,763 

4* 

35.00 

100.00 

31 ,725 

1 5 .05 

124.57 

20,766 

3 

45.00 

160.00 

37,188 

15.91 

157.97 

22,064 

2 

40.00 

140.00 

35,303 

20.97 

156.06 

28,450 


* Optinrum number of beams s 4-4 with a spacing of 376.4 cms. 



saadnj pyusnou^ ui isoo uunuji 
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freedom. The spacing of beams is 500 cms. Hence the slab 
thickness turns out to be 12.5 cms, the minimum to be 

O 

provided . Thus the dead weight or, the grid is 300 Kg/m'^ . 

All other details are worked out os explained in example 1 . 
This problem turns out to have 2 'esign variables and 39 
constraints . 

Taking the starting point as (40.0, I 4 O.O) having 
a totall- cost of Es. 1 ,33,396 the proposed optimum design 
obtained 3.s having a total cost of Rs, 79,745. 

A parametric study with respect to number of 
beams is made in the range Fb equal to 5 to 12. Table 2 
shows the results of this parametric study. It is seen from 
Table 2 that the optimum number of beams is 6 having a 
spacing between them equal to 444.44 cms. The total cost 
of roof af the optimum number of beams is Rs. 73,340. The 
parametric study of the optimal design taking Hb as a para- 
meter j.s represented graphically in Figure 5.2. 

5.4 ILLHSTRATIVE EXAECPEE 3 

-4 third grid with overall size 40 H x 40 M (Grid 3) 
is studied for optimization with respect to the spacing 
of beams. This grid is optimized for the number of beams 
Fb in the range of 5 to 12. The details for each case are 
worked out as explained in previous examples. The results 
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Table 2 



Grid 2 

Size 30 M X 30 IT 


Fo . of 

beams in each 
direction Fb 

? 

Starting design point | 

_ t 

Proposed optimum design point 

b cms 

! 

1 

}D cms 
! 
t 

T 

t 

Cost of J 
roof J 

(Rs. ) ; 

f. - .! 

b cms 

! 

D cms ! 

T 

! 

Cost of 
roof 
(Rs. ) 

9 

45.00 

200.00 

1 ,74,928 

15.01 

172.47 

79,905 

8 

50.00 

230.00 

1 ,86,093 

15.00 

189.78 

76,495 

7 

50.00 

225.00 

1 ,69,982 

15.07 

224.25 

73,864 

6* 

50.00 

225.00 

1 ,59,510 

15.04 

225.12 

^3,346 

5 

40.00 

140.00 

1,33,396 

16.67 

246 .80 

79,745 

4 

55.00 

260.00 

1,64,356 

20.00 

259.00 

90,961 

3 

55.00 

260.00 

1,56,864 

25.12 

320.56 

1 ,07,068 

* Optimum number of 

beams ; 

6-6 with a spacing of 428.57 cms. 




Optimum cost 
of roof 


Optimum number 
of beams = 6 


Optimum cos 
of grid 




RiATION IN THE OPTIMUM COST OF ROC 
ST OF GRID WlTHiRESPECT TO JBE CH 
BEAMS IN EACH DIRECTION FOR GRID 
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of the paxametric study with respect to number of beams are 
given in Table 3. The variation of optimum cost with the 
mumber of beams in each direction is shown graphically in 
Figure 5.3* 

From Table 3 it is seen that the optimum number 
of beams Hb leading to the minimum cost design is 8 at 
which the beams are at a spacing of 4-44.44 cms. 


5.5 EBSUITS 

Is already stated the total cost of roof at the 
starting point and at the imposed optimum design for 
different values of Mb are given in Tables 1 , 2 and 3 for 
G-rids1 5 2 and 3 respectively. Also the variation in the 
optimum cost with Fb for Grids 1 , 2 and 3 is shown graphi- 
cally in Figures 5.1, 5.2 and 5.3 respectively. 

The detail results of the study of optimal designs 
for Grids 1, 2 and 3 are presented in Tables 4, 5 and 6 
respectively. Table 4 gives the results for Grid 1 showing 
the variation in the optimum cost of roof, optimum cost of 
grid and the cost of slab as the number of beams vary 
between 2 and 9, Tables 5 and 6 give similar results for 
Grids 2 and 3 when lb varies in the range of 3 to 9 and 
5 to 12 respectively. 


CENi.s-... : 


5iU$ 
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Table 3 



arid 3 ; 

Size 40 1 

I X 40 M 



Fo . of 

jstarting design point 

1 - - 

i 

ilroposed optiiiium design point 

beams in each 
direction Fb 

b cms 

T 

|D cms 

1 

f 

f 

-1 

Cost of 
roof 
(Rs.) 

\ 1 

1 1 

■ b cms 1 

f 

D cms 

1 Cost of 
! roof 
j (Rs. ) 

12 

50.00 

320.00 

5,00,142 

15.06 

327.55 

1 ,95,348 

11 

45.00 

190.00 

3,91 ,109 

15.06 

320.49 

1,85,512 

10 

50.00 

310.00 

4,42,023 

15.01 

318.75 

1,76,179 

9 

55.00 

360.00 

4,99,841 

15.09 

336.00 

1,73,294 

8* 

60.00 

350.00 

4,65,450 

15.09 

348,38 

1,70,441 


60.00 

350.00 

4,46,593 

16.73 

349.22 

1,80,655 

6 

65.00 

335.00 

4,46,118 

19.19 

341 .60 

1,95,031 

5 

60.00 

340.00 

4,00,447 

22.30 

359.64 

2,18,732 

* Optimum number of 

beams; 8-8 with i 

a spacing 

of 444. 

44 cms . 




Details at the optimum design points ; Grid 1 j size 18.82 M x 18,82 M 
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Details at the optiinum design points ; Grid 
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Tables 4, 5 and 6 filso giye the weight carried by 
the grid per "unit area and '.'he macimum deflection in the 
grid for optimum design points at each spacing studied for 
Grids 1, 2 and 3 respectivej.y . 

The design details at th^^ optimum sections giving 
the average percentage of tension and compression reinfor- 
cements as well as shear and torsional steel are presented 
in Tables 1 , 2 and 3 . 

5.6 DISCUSSION AND CONCLUSIONS 

It is seen from Tables 3> 4 and 5 that for the 
spacing of beams less than or equal to 375 cms the slab 
cost is constant for a given overall size of grid. This is 
because at all the spacings less than or equal to 375 cms 
the slab thickness is 10 cms and it is provided with minimum 
specified reinforcement. Thus for 1 < 375 cms the vari- 
ation in the total minimum cost is entirely dependent upon 
the variation in the minimura cost of grid. This can be 
observed from Digures 5.1? 5.2 and 5.3 which show a constant 
difference in the ordinates of the minimum, cost of roof and 
the minimum cost of grid for spacings less than 375 cms. 

Nor a spacing greater than 375 cms upto and including 
400 cms the slab thickness remains same but the steel 
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requirement is more than the minimum specified steel. 

Hence the cost of slah increases slightly. This can be 
obserTed from the costs of slab at Hb equal to 5 and 4 for 
G-rid 1 in Table 4 and at Hb equal co 10 and 9 for Grid 3 
in Table 6. Thus the dead weight on grid is constant for 
1 :^400 cms . 

When the spacing of beams is greater than 400 cms 
the slab thickness increases which naturally causes 
increase in the cost of slat as well as in the dead load 
on grid. The slab cost increases rapidly as the spacing 
increases. The variation in the slab thickness and the cost 
of slab as the spacing becomes greater than 400 cms can be 
seen in Tables 3» 4 and 5. figures 5.1 > 5.2 and 5.3 also 
show the wariation in the cost of slab for 1^400 cms 
from the difference in the ordinates of minimum total cost 
of roof and minimum cost of grid. Thus the variation in 
the minimum total cost of roof for lr=^400 cms is also 
dependent upon the change in the cost of slab along with 
the change in the minimum cost of grid . 

To study the variation in the cost of grid as the 
number of beams in each direction changes let us see 
figures 5.1, 5.2 and 5.3. The figures show that for 
1^475 cms the cost of grid increases as the spacing 
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decreases and for 1 >475 cms the cost of grid increases as 
the spacing increases. The requirement of the minimum width 
of beam increases as 1 is increased beyond 450 cms and for 
1 ■===:450 cms the minimum width of beam is constant. It is 
seen from the optimum design points at various spacings that 
minimum specified width is always reached at the optimum. 
Thus the larger width requirement for 450 cms gives more 
dead weight of the grid which otherwise would have been 
reduced. Moreover as already explained the dead weight of 
slab also increases rapidly as the spacing increases beyond 
400 cms. Thus the total load carried by grid at larger 
spacings is more. This is observed in Tables 4, 5 and 6 
where the load per unit area carried by grid at the optimum 
design points is given for Grids 1 , 2 and 3 respectively. 

Tables 4, 5 and 6 also show that at the optimum 
designs for l'<445 cms the total load carried by the grid 
increases as the spacing decreases. This is obviously due 
to the increase in the total number of beams in the grid 
as the weight of the slab remains constant at all these 
spacings . 

From this it can be concluded that the optimum 
beam spacing leading to the minimum cost design also turns 
out to be a minimum weight design. 
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Regarding the sections of heams at the optiraum 
design for each spacing it can be said that always the 
miniimim, specified width is obtained . Though the variation 
in the depth of sections does not show any fixed trend 
still it can be observed that at larger spacings relatively 
more depth is obtained. 

looking at Tables 3? 4 and 5 the variation in the 
percentage of tension steel shows that at larger spacings 
the percentage of tension steel is higher. At the spacing 
decreases the percentage of tension steel also reduces. 

It was observed from the detail design of the optimum 
points that at lower spacings only the sections nearer to 
the supports and those near the centre of the grid require 
more steel and all other sections are required to be 
provided with the minimum specified reinforcement. Prom 
the torsional reinforcements for optimal designs given in 
Tables 4, 5 and 6 it is seen that the geometry of the 
section at the optimum requires very small amount of 
torsional steel. 

The maximum deflections at optimum designs given 
in Tables 4, 5 and 6 show that the order of maximum defle- 
ction at optimum design points is very small. 

Prom the above discussion it can be concluded 
that the optimum spacing of beams is the one which is 
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aroimd 400 cms for grids up to a span of 40 M. The minimum 
cost design at the optimum spacing also turns out to be a 
minimum weight design. It is also seen that at optimum 
designs the torsional steel requirement is very small, 
further the optimal design gets to the minimum specified 
width of the beam and the deflections at the optimal designs 
are of small order. 

Prom the optimum number of beams obtained for the 
three grids i.e. Grids 1, 2 and 3 } an empirical formula 
for obtaining the optimum number of beams in a grid of 
given size and of the type considered in the present 
problem can be suggested. The proposed formula is as 
follows 

jr. _ Span of grid in meters 

^optimum “ 5 

Prom the above formula the number of beams at optimum 
design of Grids 1 , 2 and 3 respectively turn out to be 4, 

6 and 8 as seen in Tables 1 , 2 and 3 respectively. 
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